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ABSTRACT 

We present an analysis of the atmospheric content of aerosols measured at Observato- 

rio del Roque de los Muchachos (ORM; Canary Islands). Using a laser diode particle 

counter located at the Telescopio Nazionale Galileo (TNG) we have detected particles 

of 0.3, 0.5, 1.0, 3.0, 5.0 and 10.0 /im size. 

The seasonal behavior of the dust content in the atmosphere is calculated. The Spring 

has been found to be dustier than the Summer, but dusty conditions may also occur 

in Winter. 

A method to estimate the contribution of the aerosols emissivity to the sky brightness 

in the near-infrared (NIR) is presented. The contribution of dust emission to the sky 

background in the NIR has been found to be negligible comparable to the airglow, 

with a maximum contribution of about 8-10% in the K$ band in the dusty days. 

Key words: Site testing - Methods: data analysis - Methods: statistical. 



1 INTRODUCTION 

Superb observing conditions are crucial in order to obtain 
the best scientific output form the next generation of ground 
based telescopes, and this requires monitoring all relevant 
parameters that may affect observations. Due to their po- 
tentially detrimental impact on astronomical observations, 
atmospheric aerosols are among the most important param- 
eters to be monitored in modern site testing campaigns. The 
performance and the safety of telescopes depend on the pres- 
ence of atmospheric dust, which may deposit on mirrors, 
increase atmospheric extinction, and emits in the infrared 
(IR) bands thus increasing the sky brightness. 

Several studies about the atmospheric radiative effects 
of mineral dust exist. Gelado et al. (2003) analysed the mean 
aerosols content at Gran Canaria (Canary Islands) in peri- 
ods 1997-1998 and 2002-2003. A mean grain size of 0.6-4.9 
/im and a large annual variability in both density and size 
distribution have been found. 

Cuevas&Baldasano (2009) provided information on the 
incidence of the dust loaded African air masses at Observa- 
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torio del Teide in Tenerife (Canary Islands) and at Obser- 
vatorio del Roque de Los Muchachos (ORM) in La Palma 
(Canary Islands) . The analysis has been made using the high 
quality observations performed by the Izafia atmospheric ob- 
servatory. The in situ measurements for the 2002-2008 pe- 
riod show that the air is only partially affected by some dust 
loaded by African air mass intrusion. A significant PM10 
concentration has been found only above the 80 percentile 
(26 fig m -3 ) at Izafia. 

The first analyses of the optical properties of d ust an d 
their impact on astro nomy were done by M urdin (|1985T) . 
Stickla nd et al. |l987j), Guerrero et al. (1998J), and J imenez 
et al. (119981 ). More recently, Lombardi et al. (120081 ) (here- 
after Paper I) analyzed 5+ years of time-series data span- 
ning the period August 2001 to December 2006 obtained 
with the Telescopio Nazionale Galileo (TNG) dust counter, 
and showed that dust particles increase the extinction in the 
B, V and / bands at ORM. 

In the present paper we extend the analysis of Paper 
I by adding almost 3 years of new data (March 2007 to 
January 2010) to cover a total period of almost 8 years. In 
the first part of the paper we calculate the seasonal trend 
of the atmospheric aerosol content on yearly and monthly 
basis. In the second part of the paper we present a method 
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Table 1. Main characteristics of the dust monitors. 



ABACUS TM301 



LASAIRII 310B 



input flow rate 

size channels 

light source 

sample rate 

output 

running time 



0.1 Cubic Foot Minute 

0.3, 0.5, 1.0, 5.0 fim 

laser diode (A = 780 nm) 

1 data per minute 

RS-232 

August 2001 - December 2006 



1.0 Cubic Foot Minute 

0.3, 0.5, 1.0, 3.0, 5.0, 10.0 /im 

laser diode (A = 780 nm) 

1 data per 2 hours 

Ethernet 

March 2007 — today (January 2010 for this paper) 



to estimate the thermal emissivity of the dust in the near- 
infrared (NIR). 



2 DUST MONITORS 

Since 2001 the TNG site monitoring group has used two 
different particle counters made by Particle Measuring Sys- 
tem to monitor the dust content of the atmosphere around 
the telescope. The Abacus TM301 measured between Au- 
gust 2001 and December 2006, and the Lasairll 310B is in 
operation since March 2007. 

The data from the Abacus TM301 were extensively ana- 
lyzed and presented in Paper I. The Lasairll 310B provides a 
significant improvement over its predecessor, in fact it mea- 
sures 6 different particle sizes (0.3, 0.5, 1.0, 3.0, 5.0 and 10.0 
/im) instead of only 4 measured by the Abacus TM301. 

Table Q] summarizes the basic instrumental characteris- 
tics of these counters that use a laser scattering technique 
for environmental ambient air analysis. Both are compact 
and portable devices designed to measure the air-purity 
of closed envi ronm ents such as clean rooms (Porce ddu et 
al.(author?) |2002| ; Ghedina et al. (author?) 12004 Paper 
I). The light scattered by the dust particles is converted 
to voltage pulses of amplitude proportional to particle size 
and frequency proportional to particle density. A long silicon 
pipe through the enclosure feeds the pump with external air 
at 13 m above the ground, that correspond to the level of 
the TNG primary mirror. 



3 SEASONAL DISTRIBUTION OF THE DUST 
IN BACKGROUND AND DUSTY 
CONDITIONS 

The raw data from the counters have been analyzed follow- 
ing the same method of Paper I. Because we cannot distin- 
guish between the types of particles, but only their size, we 
rejected data obtained when the measured relative humid- 
ity was greater than 85%, that corresponds to the conden- 
sation point for water vapor particles. As we will see later, 
this may introduce biases in the most humid months. Dust 
counts having a values of a few a above the median values 
have been classified as dust storms. The background is eval- 
uated using a k — a c lipping algorithm as described in Huber 
(|l98ll) ; Patat (|200St ); and also used in Paper I. 

Because particle counts Ni -where i is the ith size bin- 
have a Poissonian distribution, for each fixed period (year, 
semester, or month) we compute the median (M), and for 
each particle size we calculate the Median Absolute Devia- 
tion (P) defined as the median of the distribution \Ni — M\. 



For a Gaussian distribution, the ratio between the standard 
deviation and T is equal to 1.48 (Huber(author?) 1981), so 
we set a = 1.48P. We perform two iterations to compute the 
background rejecting in each iteration counts that exceed, 



\Ni - M\ > K<7 



(1) 



In the first iteration we set k — 3 and in the second we 
rejected measurements between —3a and +2a in order to 
finally get the background distribution, uncontaminated by 
dust peaks corresponding to dust storms. 

In order to investigate the presence of seasonal vari- 
ations of the background dust content, we divided our 
database in wintertime (semester October-March) and sum- 
mertime (semester April-September). Figure Q] shows the 
trend of the annual background distribution for the 8- years 
measured particles (0.3, 0.5, 1.0 and 5.0 /j,m). We remark 
that in this plot we have removed dust storms using the 
K — a clipping procedure described above, so these results 
correspond to strictly background dust conditions. Figure [T] 
shows that the dust background is always larger in Summer 
than in Winter except in 2002, when the largest particles 
(5.0 /im) were significantly more abundant in winter than 
in summer. Curiously, this anomaly is not present for the 
smaller particles. A further comparison with the TOMS 
(Total Ozone Mapping Spectrometer) data does not show 
the same result. It is important to mention that the TOMS 
aerosol index (AI) has as a major limitation the inability 
to detect dust occurring at or near the surface since the 
ground signal can o verwhelm the dust signal (Her man e t 
al. (author?) 1997; Gao&Washington(author?) 2010T I. 
The 2002 dust recorded by our ground based dust sensor 
could be confined to the lowermost layers of the atmosphere, 
which might explain why it could not be seen by TOMS. 
Our guess is that the 2002 phenomenon corresponds to a 
local event not seen by the satellite (contamination from 
local recycled dust?). 

Table [2] presents the dust background in counts per cu- 
bic meter [N m -3 ] in background dust conditions distin- 
guishing between wintertime, summertime, and an entire 
annual cycle, and also the dust content during dust storms 
in an entire annual cycle, for 3 different cases: 

(1) using the Abacus TM301 database between August 
2001 and December 2006 for particles of 0.3, 0.5, 1.0, and 
5.0 (im sizes; 

(2) using the Lasairll 310B database between March 2007 
and January 2010 for particles of 0.3, 0.5, 1.0, 3.0, 5.0 and 
10.0 /an sizes; 
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Table 2. Dust content in [N m 
text in Section l3t . 



at ORM in wintertime, summertime, and in the entire annual cycle for cases (1), (2) and (3) (see 



CASE 


(1) 




BACKGROUND 




DUSTY 




statistics 






Wintertime 


Summertime 


Annual 


Annual 




running time 


0.3 


/jin 


1.3 X 10 6 


4.4 X 10 6 


3.0 x 10 6 


1.8 X 10 7 


Aug. 


2001 - Dec. 2006 


0.5 


/Jin 


1.2 X 10 5 


3.8 X 10 s 


2.5 x 10 5 


6.1 x 10 6 




5J 


1.0 


jLtm 


0.5 x 10 5 


1.5 x 10 5 


1.0 x 10 5 


4.1 x 10 6 




" 


5.0 


/tm 


0.7 x 10 3 


1.5 x 10 3 


1.1 x 10 3 


1.6 x 10 5 




" 




CASE 


(2) 




BACKGROUND 




DUSTY 




statistics 






Wintertime 


Summertime 


Annual 


Annual 




running time 


0.3 


/tin 


1.1 X 10 6 


3.7 X 10 6 


2.3 x 10 6 


3.0 x 10 7 


Mar 


2007 - Jan. 2010 


0.5 


/Jin 


1.2 X 10 5 


3.7 x 10 s 


2.5 x 10 5 


4.1 x 10 6 




5J 


1.0 


/(in 


0.5 x 10 5 


1.0 x 10 s 


0.5 x 10 5 


2.6 x 10 6 




» 


3.0 


/(in 


2.7 x 10 3 


5.2 x 10 3 


3.9 x 10 3 


5.1 x 10 5 




" 


5.0 


/(in 


1.2 x 10 3 


1.6 x 10 3 


1.5 x 10 3 


1.6 x 10 5 




" 


10.0 


/(in 


1.2 x 10 2 


0.7 x 10 2 


1.0 x 10 2 


1.1 x 10 4 




" 




CASE 


(3) 




BACKGROUND 




DUSTY 




statistics 






Wintertime 


Summertime 


Annual 


Annual 




running time 


0.3 


/(in 


1.0 X 10 6 


3.7 X 10 6 


2.3 x 10 6 


2.8 X 10 7 


Aug 


2001 - Jan. 2010 


0.5 


/Jin 


1.2 X 10 5 


3.7 x 10 5 


2.4 X 10 5 


4.6 x 10 6 




» 


1.0 


/jm 


0.4 X 10 5 


1.1 x 10 s 


0.6 x 10 5 


2.6 x 10 6 




5J 


3.0 


/jm 


2.7 x 10 3 


5.2 x 10 3 


3.9 x 10 3 


5.1 x 10 5 


Mar 


2007 - Jan. 2010 


5.0 


/Jin 


1.0 x 10 3 


1.5 x 10 3 


1.3 x 10 3 


1.6 x 10 5 


Aug 


2001 - Jan. 2010 


10.0 


/(ill 


1.2 x 10 2 


0.7 x 10 2 


1.0 x 10 2 


1.1 x 10 4 


Mar 


2007 - Jan. 2010 
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Figure 1. Seasonal dust background distribution at TNG during background dust conditions, in counts per cubic meter [N m 3 ]. 
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(3) using both the Abacus TM301 and the Lasairll 310B 
databases between August 2001 and January 2010 for parti- 
cles of 0.3, 0.5, 1.0 and 5.0 /im sizes, and the Lasairll 310B 
database for particles of 3.0 and 10.0 /xm sizes. 

Figure [JJ and Table [2] show that the prevailing aerosol 
conditions at ORM have not changed significantly during 
the 8 years of monitoring, with summertime being dustier 
than wintertime, with the exception of the 10.0 (j,m particles 
that show the opposite behavior. 

We now consider the monthly median counts without 
excluding the dust storms. In Paper I we found that be- 
tween 2001 and 2006 the monthly distribution of aerosols 
is characterized by a significant increase in the counts dur- 
ing February-April and July-August of each year. Figure 
2 shows the median monthly distribution of particles from 
the Lasairll 310B data that confirms the trends found in 
Paper I, with the exception of November 2007 that was 
clearly a very dusty month. June, September, October, and 
December- January appear to be the cleanest months, but 
we remark that December and January are also the most 
humid months, so these results may be biased because the 
dust counters are not reliable (relative humidity greater than 
85%). 

Figure [3] shows the monthly distribution of overall dust 
density Mio (in micro grams per cubic meter) computed 
adding together all particle sizes for Case (2) and assuming 
that all particles have the same density of 2.5 g cm -3 , which 
is typ ical of silicates and quartz aggregates (Suh(author?) 
Il999f ) that are the m ain components of Saharan dust (Mur- 
din(author?) 1985). In this plot July emerges as the dusti- 
est month, while somewhat surprisingly June is the cleanest. 
As already discussed, the Spring is dustier than the Summer 
and dusty conditions are also frequent in the Winter. 



4 THERMAL BACKGROUND IN THE NIR 

A dust particle of radius a is heated by an ambient gas of 
molecules at temperature T m at a rate that depends on the 
number density of molecules and the fraction of energy E 
that is deposited i n the dust grain by the impinging particles 
(Dwek(author?) Il986h . The dust particle radiates in the 
infrared the energy acquired in the collision at a rate that 
depends on the grain temperature Td as, 



E* = Ana 



nB x (T d )Q x (a)d\ 



(2) 



where B\(Td) is the Planck function for temperature Td and 
Qx (a) is the grain emissivity at wavelength A that i s usually 
calculated using Mie theory (Dewk(author?)|l986|). In this 
Section we estimate the thermal background emission of at- 
mospheric dust in the NIR spectral bands, for which we will 
use the 2MASS photometric system reproduced in Table [3] 
We can use our measured dust densities to calculate the 
total dust optical depth t\ as the sum of the contribution 
of the optical depth of each particle size, 



r\ = ^V A (a) = 






N(a)a x {a) 



(3) 



where N(a) is the column density of particles of radius a, and 
a\(a) is the wavelength dependent absorption cross section 
from the Mie theory for particles of radius a. 
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Figure 2. Seasonal distribution of the monthly dust at ORM in 
counts per cubic meter [N m — 3 ], Different years are represented 
by different symbols. The solid line shows the monthly values 
averaging over all years. 
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Figure 3. Monthly distribution of dust density at ORM com- 
puted adding together all particle sizes for Case (2). A constant 
density of 2.5 g cm -3 was assumed for all particle sizes. 



Table 3. Absolute calibration of the 2MASS photometric system 
taken from Cohen et al. {.2003). 



Filter ^eff 
[H 



AA 

H 



[W 



n-2 



/im 



J 
H 
K S 



1.235 
1.662 
2.159 



0.162 
0.251 
0.262 



3.129 x 10 -9 
1.133 x 10 -9 
4.283 x 10 



in 



Satellite measurements indicate that the dust density 
over the Canary Islands is approximately constant at alti- 
tudes between 2500 m and 5000 m above sea level, and drops 
to vi rtually zero above 5000 m (Smirnov et al. (author?) 
2002J : Hsu et al. (author?) 1 19991 ; Alpert et al. (author?) 
2004 ; Paper I). Thus, the column densities are simply the 
volume densities multiplied by 2500 m. The resulting optical 
depths are calculated using the dust content in background 
and dusty conditions obtained for an entire annual cycle and 
shown in Case (3) of Table [2] Results are given in Table [4] 
and are in very good agreement with those in the visible from 
Smirnov et al. (|2002r i calculated above 2356 meters above sea 
level in Tenerife Canarian Island (see Table 2 in the men- 
tioned paper). The extrapolation to NIR wavelengths still 
remains in good agreement with our results. 



Table 4. Aerosol optical depth calculated for the 2MASS filters. 



Filter background dusty 



J 
H 
K S 



0.006 
0.005 
0.004 



0.296 
0.280 
0.270 



Table 5. Aerosols background NIR emissions in [mag arcsec 2 ] 
in background and dusty days. 

Filter background conditions dusty days 



./ 
H 
K S 



39.2 
27.8 
20.3 



35.2 
23.6 
15.8 



The cooling times for dust particles of the sizes and 
temperatures typical of the Saharan dust above OR M are 
very short (Kais er(au thor?) Il970l ; Draine(author?) Il98ll ; 
Dewk(author?)ll98y), so it is reasonable to assume that the 
dust is locally in thermal equilibrium with the surrounding 
air, while clearly the dust is optically thin even under the 
dustiest conditions. Thus, the thermal radiation emitted by 
the dust particles can be calculated as, 



h=T X B x (T d ) = 



2c 2 /i. 



A 5 exp( 7 



(4) 



where B\(Td) from Tokunaga (|2000l ) is expressed in [W m -2 
pan -1 sr -1 ], c is the light speed in vacuum, h is the Planck's 
constant, and ks is the Boltzmann constant. 

We must calculate the Planck function considering that 
Td decreases with altitude (h) from the median value at 
the g round (about 282 K, see Lombardi et al. (author?) 
120061 ) with a wet vertical adiabat ic lap se rate of —0.006 K 
m" 1 (Kittel&Kroemer(author?) 1980). Thus, the thermal 
spectrum of the dust is given by the expression, 



F x 



r 2 

Jo 



B x (T d (h))dh 



(5) 



where the sum is over all particle sizes, a. Using the absolute 
calibration of the 2MASS photometric system presented in 
Table [3] we obtain the aerosols sky brightness in the NIR in 
background and dusty days as, 



m\ = —2.5 log 



Fa 



(6) 



The results are shown in Table [5] that clearly show that the 
sky brightness due to dust in the atmosphere is significant 
only during dusty conditions, and only in the Ks band. From 
the TNG data archive we know that the sky background in 
J is between 15.0 and 16.0 mag arcsec -2 , in H is between 
13.4 and 14.7 mag arcsec -2 , while in Ks it varies from 12.5 
and 13.0 mag arcsec -2 . We conclude that at ORM the con- 
tribution of the dust to the sky background in the NIR is 
mostly negligible in both background and dusty conditions, 
with a maximum contribution of about 8-10% in the Ks 
band in the dusty days. 

Given that our results were computed over a long pe- 
riod, the averaging process could have smoothed larger 
events, so the full radiative impact of these events may go 
undetected. As an example, we have calculated the optical 
depth (and thus the resulting thermal emissivity) for a ma- 
jor event occurred between 25 and 30 July 2007. Results are 
reported in Table |6] and confirm that still during a major 
event the contribution of the dust to the sky background in 
the NIR is not critical. 
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Table 6. Aerosols optical depth and NIR emissions in [mag 
arcsec -2 ] for the major event occurred between 25 and 30 July 
2007. 



Filter 


T\ 


emission 


J 


0.388 


34.5 


11 


0.372 


23.0 


K S 


0.356 


15.2 



We can conclude that the obtained thermal emission is 
the contribution of our in situ measured particles and results 
could be intended as a lower limit. 



5 CONCLUSIONS 

Using a laser diode particle counter near the TNG at ORM 
we have measured the densities of airborne aerosols of 0.3, 
0.5, 1.0, 3.0, 5.0 and 10.0 /an size. The seasonal trends of 
the particle content in the atmosphere have not changed sig- 
nificantly between 2001 and 2010. The monthly distribution 
of aerosols is characterized by an increase during February- 
April and July-August of each year: the Spring is dustier 
than the Summer, but dusty conditions may also occur in 
Winter. 

Using the Mie theory we calculated the dust absorption 
cross section and thus estimated the thermal emission in the 
2MASS NIR spectral bands. Assuming that the dust parti- 
cles are locally in thermal equilibrium with the surrounding 
air, and that the air temperature decreases with altitude 
with the wet vertical adiabatic lapse (—0.006 K m _1 ), we 
found that the contribution of dust emission to the total 
sky background in the NIR is negligible comparable to the 
airglow component during both background and dusty con- 
ditions, with a maximum contribution of about 8-10% in the 
Ks band in the dusty days. 



Gelado, M. D., Dorta, P. J., Torres, M.E., Hernandez, J. J., 

Collado, C, Siruela, V. F., Cardona, P., & Rodriguez, M. 

J. 2003, Caracterizacion del aerosol Sahariano en Gran Ca- 

naria, Encuentro sobre Metcorologia y Atmosfera de Canarias, 

http://acceda.ulpgc.eS/bitstream/10553/1531/l/5391.pdf 

Ghedina, A., Pedani, M., Guerra, J. C, Zitelli, V., & Porceddu, 

I. 2004, SPIE, 5489, 227 

Guerrero, M. A., Garcia-Lopez, R. J., Corradi, R. L. M., 

Jimenez, A., Fuensalida, J. J., et al., 1998, New Astron. Rev., 

42, 529 

Jimenez, A., Gonzalez Jorge, H., & Rabello-Soares, M. C. 1998, 

A&AS, 129, 413 

Herman, J., Bhartia, P., Torres, O., Hsu, C, Seftor, C, &; Celar- 

ier, E. 1997, J. Geophys. Res., 102(D14), 16911 

Huber, P. J. 1981, Robust Statistics, ed. J. Wiley & Sons, NY 

Hsu, N. C, Herman, J. R., Torres, O., Holben, B. N., Tanre, D., 

& Eck, T. F. 1999, J. Geophys. Res., 104, 6269 

Kaiser, C. B. 1970, ApJ, 159, 77 

Kittel, C, & Kroemer, H. 1980, Thermal Physics, ed. W. H. 

Freeman Company (2nd ed.) 

Lombardi, G., Zitelli, V., Ortolani, S., & Pedani, M. 2006, PASP, 

118, 1198 

Lombardi, G, Zitelli, V., Ortolani, S., Pedani, M., & Ghedina, 

A. 2008, A&A, 483, 651 (Paper I) 

Murdin, P. 1985, Vistas Astron., 28, 449 

Patat, F. 2003, A&A, 401, 797 

Porceddu, I., Zitelli, V., Buffa, F., & Ghedina, A. 2002, SPIE, 

4844, 358 

Smirnov, A., Holben, B. N., Slutsker, I. , Welton, E. J., & For- 

menti, P. 1998, J. Geophys. Res., 103(D21), 28, 079-28, 092 

Suh, K. W. 1999, MNRAS, 304, 389 

Stickland, D. J., Lloyd, C, Pike., C. D., & Walker, E. N. 1987, 

The Observatory, 107, 74 

Tokunaga, A. T. 2000, in Allen's Astrophysical Quantities, ed. 

A. N. Cox (4th ed.; New York: AIP) 

This paper has been typeset from a T^X/ P/IpX. file prepared 
by the author. 



ACKNOWLEDGMENTS 

The authors aknowledge the anonymous reviewer for the 
very helpful comments on the paper. 



REFERENCES 

Alpert, P., Kishcha, P., Shtivelman, A., Krichak, S. O., & 
Joseph, J. H. 2004, Atmos. Res., 70, 109 

Cohen, M., Wheaton, W. A., & Megeath, S. T. 2003, AJ, 126, 
1090 

Cuevas E., & Baldasano J. M. 2009, Report on the Incidence of 
African dust intrusions at the Astronomical Observatories ofteh 
Canary Islands: characterization and temporal analysis, Report 
released by the Izana Atmospheric Research Center (Meteoro- 
logical State Agency of Spain; AEMET) and the Earth Sciences 
Department (Barcelona Supercomputing Center - Centro Na- 
tional de Supercomputacion; BSC-CNS) under request of the 
Instituto de Astrofisica de Canarias (IAC) 
Dwek, E. 1986, ApJ, 302, 363 
Draine, B. T. 1981, ApJ, 245, 880 

Gao, H., & Washington, R. 2010, Climate Dynamics, 35(2-3), 
511 



